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SOME  CONSIDERATIONS  ON  AN  AIRFOIL  IN 
AN  OSCILLATING  STREAM 
By  J.  Mayo  Greenberg  , 


.  SUMMARY 


The  velocity  potential,  lift  force,  moment,  and  propulsive 
force  on  a  two-dimensional  airfoil  in  a  stream  of  periodically 
varying  angle  of  attack  have  "been  derived  on  the  "basis  of  non¬ 
stationary  incompressible  potential -flow  theory  which  includes 
the  effect  of  the  continuous  sheet  of  vortices  shed  from  the 
trailing  edge .  Application  of  these  results  was  made  in  an, 
analysis  of  the  variation  with  frequency  of  the  propulsive  force 
on  an  airfoil  in  an  oscillating  stream  and  in  an  analysis  of  the 
problem  of  forced  vibrations  of  an  airfoil  in  an  oscillating 
stream  with  consideration  of  the  stiffness  of  the  airfoil  and 
the  position  Qf  its  torsion  axis..  It  was  shown  that  when  the 
torsion  axis  of  the  airfoil  is  ahead,  of  the  quarter- chord  point 
the  amplitude  of  vibrations  is  generally  not  large,  but  when  the 
torsion  axis  is  behind  the  quarter-chord  point  certain  conditions 
exist  under  which  dangerous  amplitudes  of  vibration  may  occur. 

The  nonuniform  response  which  was  found  for  a  freely  hinged 
airfoil  restricts  the  use  of  such  a  device  as  a  flow-measuring 
Instrument  to  the  measurement  of  only  very  low-frequency  angular 
variations  in  an  oscillating  stream. 

It  is  expected  that  the  results  of  the  theoretical  treatment 
of  the  propulsive  force  will  be  useful  in  considerations  of 
counterrotating-propeller  efficiencies  and  that  the  analysis  of 
the  problem  of  forced  vibrations  will  be  useful  in  design  con¬ 
siderations  of  lifting  surfaces  operating  in  oscillating  streams; 
for  example ,  wind-tunnel  fan  blades  behind  a  set  of  prerqtation 
vanes,  'or  tail  surfaces  in  fluctuating  wakes. 


INTRODUCTION 


The  phenomenon  of  an  airfoil  in  an  oscillating  stream  (that 
is,  a  stream  of  which  the  angle  of  attack  varies  periodically)  is 
encountered  in  many  phases  of  aeronautics.  For  example,  the 
effect  of  a  set  of  prerotation  vanes  upon  a  wind-tunnel  fan  blade 
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is  to  produce  periodic  disturbances  through  which  the  fan  blades 
pass.  Further,  a  consideration  of  counterrotating  propellers 
shows  that  the  rearward  blades  operating  in  the  helical  wake 
produced  by  the  forward  blades  are  in  a  stream  of  varying  angle 
of  attack.  As  another  example,  the  horizontal  tail  of  an  air¬ 
plane  may  be  subjected  to  forces  induced  by  fluctuations  in  the 
angle  of  the  wing  wake .  • 

The  theory  of  nonstationary  motion  around  airfoils  with 
consideration  of  partial  motions  of  the  fluid  has  been  developed 
by  Garrick  (reference  1),  Kiissner  (reference  2),  and  others.  The 
present  alternative  treatment  leads  to  a  value  for  the  lift  force 
on  an  airfoil  in  an  oscillating  stream  which  is  in  agreement 
with  that  obtained  in  references  1  and  2.  The  treatment  in 
reference  2  leads  to  general  results  for  the  propulsive  force. 

The  present  paper,  employing  somewhat  different  derivations, 
gives  explicit  results  for  the  propulsion  as  well  aB  for  the 
lift  on  an  airfoil  in  an  osoillating  stream. 

Two  problems  which  arise  In  cases  of  oscillating  flows 
are  (l)  the  production  of  vibrations  and  (2)  the  so-called 
"Katzmayr  effect"  (reference  3)  or  existence  of  a  propulsive 
force.  Vith  regard  to  problem  (l)  the  object  of  the  present 
paper  is-  to  examine  theoretically  the  dynamics  of  an  airfoil  in 
an  oscillating  stream  and  to  determine  unde_r  what  conditions 
dangerous  amplitudes  of  vibration  may  occur.  A  special  case 
for  which  the  torsional  stiffness  is  zero  is  treated  with  a 
view  to  the  possibility  of  using  a  small  freely  hinged  airfoil 
as  a  device  for  measuring  the  angular  amplitude  of  an  oscillating 
stream.  Previous  work  on  this  problem  (reference  1)  has  been 
done  for  the  case  in  which  the  stiffness  of  the  airfoil  was 
expected  to  give  large  vibrations .  Vith  regard  to  problem  (2), 
which  is  of  importance  in  considerations  of  counterrotating- 
propeller  efficiencies,  a  theoretical  investigation  is  made  of 
the  horizontal  forces  experienced  by  an  airfoil  in  an  oscillating 
stream. 

The  theoretical  development  is  divided  into  three  parts: 

(1)  derivation  of  the  lift  force  and  moment  acting  on  an  airfoil 
in  an  osoillating  stream,  (2)  derivation  of  the  propulsive  force, 
and  (3)  derivation  and  solution  of  the  equation  of  motion  of  an 
airfoil  executing  torsional  vibrations  in  an  oscillating  stream. 
The  theoretical  methods  used  in  the  derivation  of  the  lift  forces 
and  moments  consist  in  an  extension  of  the  methods  of  Theodorsen 
(reference  5)*  For  the  derivation  of  the  propulsive  force 
application  is  made  of  the  method  outlined  by  von  Kerman  and 
Burgers  (reference  6  pp.  52  and  30&)  •  The  following  visual  assump¬ 
tions  are  made  throughout:  (a)  incompressible  potential  flow, 
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(b)  two-dimensional  flat-plate  airfoil,  (c)  small  oscillations, 
and  (d)  plane  wake  extending  from  trailing  edge  to  infinity. 


SYMBOLS 


a 

a 

P 

x 

t 

V 

V 

<JJ 

k 

C(k) 

F,G 

p 

P 

■P 


R« 


half  chord  of  airfoil 

x-coordinate  of  torsion  axis  of  airfoil 

angle  of  attack  of  -  airfoil '  measured  clockwise  ..from 
horizontal 

angle  of  airstream  from  mean  direction- measured 
positive  counterclockwise 

horizontal  coordinate)  nondimenaional  with  respect  to  h 
time 

,  stream  velocity 
frequency 
circular  frequency 
reduced  frequency  • 

Theodorsen’s  C-function  from  reference  5  (F  +  iG) 

Zeal  and  imaginary;-parts  of  C~function 
local  static  [pressure 

air  density  .  ■. 

perpendicular  force . 

pitching  moment  about  x  =  a  measured  positive 
counterclockwise 

moment  of  inertia,  about  x  -  a 

torsional  stiffness  of  wing 

propulsive,  force 
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<P 

% 

U 

t 


none iroulatory  velocity  potential 
circulatory  velocity  potential 
strength  of  wake  discontinuity 
phase  angle 


JQ(k) , J^fk)  Bessel  functions  of  the  first  kind  and  zero  and 
first  order 

Yo(k)  jY-^k)  Bessel  function  of,  the  second  kind  and  zero  and 
first  order 


r 

m- 

r. 


strength  of  vorticity  distribution  on  the  airfoil 
mass  per  unit  length  of  wing 
radius  of  gyration  divided  By  h 

ratio  of  mass  of  cylinder  of  air  of  diameter  equal 


to  chord  of  wing  to  mass  of  wing 
constant 


snhl 

m 


LIFT  FORCE  AND  PITCHING  MOMENT 


In  accordance  with  Theodorsen  (reference  5),  the  forces  due  to 
the  none iroulatory  flow  and  to  the  effect  of  the  wake  are  treated 
separately . 

None iroulatory  force  and  moment.-  Consider  an  airfoil  of 
chord  2Tb  at  zero  angle  of  attack  with  respect  to  the  average 
direction  of  a  sinusoidal  stream  traveling  to  the  right  (fig.  l) . 

If  the  amplitude  3q  angle-of-attack  change  in  the  stream  is 
small  then  the  horizontal  and  vertical  components  of  the  velocity, 
respectively,  are  given  by 


and 


V^  =  V  COS  P  to  v  • 


vp  =  v  sin  p  i v  vP 
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where,  with  the  assumption  of  sinusoidal  oscillations, 


3  .  p0ei^t"^ 


Thus,  the  airfoil  may  he  considered  as  being  in  a  uniform 
horizontal  stream  of  velocity  .  v  plus  a  vertical  sinusoidal 
gust  of  the  form 


e(x)  =  vPoe1^"^ 


The  velocity  potential  cp  for  such  a  normal-velocity  distribution 
is  (appendix  A) 


_  _  bYP0  (/i  .  x2-  Qi(<ot-kr) 
2c 


where  Jq(u)  is  a  Bessel  function  of  the  first  kind  and  zero 
order . 

Use  of  the  equation  of  motion  for  nonstationary  flow  gives 


Qiux  jq(u)  du 


(1) 


§w 
St  = 


where 

w  fluid  velocity 

p  local  static  pressure 

p  air  density 
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and  the  substitution  w 
at  the  point  x  as 


gives  the  pressure  difference  Ap 


ip  '  (i  1 


(2) 


Integration  of  this  local  pressure  difference  over  the  length  of 

the  airfoil  gives  a?  the •  total  force  P  (see. appendix  B,  equation  (B8)) 


P  =  gjtiphv2B0J1(k)e12>t  (3) 

where  J^(k)  is  a  Bessel  function  of  the  first  kind  ands  first 
order . 

The  none ir culatory  moment  about  x  »  a  (fig»  2}  is  obtained 
from  the  integral 


r X 

b2  /  Ap(x  -  a) 
J-l. 


dx 


which  yields  (appendix  B,  equation  (Bl6)) 


M&  =  -npb2v2p0ei<1Dt  ^2iaJ1(k)  +  JQ(k)j 


(M 


Circulatory  force  and  moment.-  The" velocity  potential  of  an 
element  of  vorticity  -AT  at  a  position  xq  in  the  wake  and  its 
mate  AF  distributed  over  the  airfoil  is  (reference  5) 
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The  element  -Ar  moves  to  the  right  relative  to  the  airfoil  with 
a  velocity  v.  Thus, 


dffxxp 

St 


v-JSBS 


SXq 


Substituting  this  expression  and 


ocp 


**0 

Sx 


into  equation  (2)  and 


UA 

integrating  the  effect  of  the  entire  wake  on  the  airfoil  yieldB  the 
force 


P  =  -pvb  /  ..^-5L,^=  TT  dsu  (6) 

ifo2  ‘  1 


where  U  cLxq  is  the  element  of  vorticity  AT  at  the  point  Xq. 

The  Kutta  oondition  requires  that  at  the  trailing  edge  of 
the  plate  the  induced  velocity  equate  zero,-  therefore,  at  x  =  1 


h  +  i 


where 


cpr  »  h  cpr^  dxQ 

Introducing  the  potential  cp  from  equation  (l)  results  in 

hf’i ’  *b  -  -  ^  ^  f*  -  (t) 
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Let' 


9  =  -  Z?2  r  eitt  J0(»)  du  (8) 

k  .  '^0 

where  u  is  an  integrating  variable.  Equation  (8)  'becomes 
(appendix  C) 


Q  *  -vp0|j0(k)  -  iJ^Ckj]  etot  (9) 

Combining  equations  (6) ,  (7),  and  (8)  and  assuming  the  wake  to  be 
of  the  form 

XT  -  xjoei(mt--laco)  .  . 

gives' for  the  circulatory  force 


F 


-2rtpvbQ 


(10) 


Similarly,  the  circulatory  moment  which  is  obtained  from 


is  (reference  5) 


Ma  =  -2rtpvbQ 


Ma  =  b2  J  Ap(x  -  a)  dx 


/ _ fo _  “ik*0 

l.fat  iWl  /*02  ~  1  8  . 

2'  \  x0  +  1  -ikxo 


dx0 


ikxo 

TVT  9 


(11) 
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The  integral  express!',  n  in  equations  (10)  and  (ll)  is  Theodorsen‘e  C- 
function  (reference  5) .  Thus, 


P  =  2rtphv%C(k)  '  J0(k)  -  iJ1(k)J  etot  (12) 


and 


Ma  =  |  -  (  &  +  £)c(k) 


J0(k)  -  Ux(k) 


.hot 


(13) 


Adding  equations  (3)  and  (12)  gives  for  the  total  force 


r 

P  ®  2irphv2(30  <jC(k)  JjQ(k)  -  iJx(k)J  +  iJ1(k) 


.hot 


(1*) 


This  expression  agrees  with  that  given  "by  Garrick  (reference  1)  and 
Kiissner  (reference  2)  in  which  somewhat  different  methods  of 
derivation  are  used  from  those  used  in  the.  present  paper. 

Adding  equations  ( 4-)  and  (13)  gives,  for  the  total  moment 
about  x  =  a. 


M 


=  -Tfp'b2v2P0  |2iaJ1(k)  +  J0(k)j 


.hot 


2irpb2v23( 


C(k) 


J0(k)  -  iJx(k) 


5hot 


(15) 


Examination  of  equations  (ll)  and  (15)  leads  to 


Ma  =  -b( a  +  1 


This  equation  means  that  the  center  of  pressure  is  at  the  quarter- 
chord  point  of  the  airfoil.. 
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■  PROPULSIVE  FORCE 

According  to  yon  Kaman-aiid  Burgers  (reference  6,  p.  52),  if 
the  strength  of  the  vortioity  distribution  at  the  leading  edge  of 
the  plate  is  of  'the  form 


then  the  suction  or  propulsive  force  acting  on  the  airfoil  a  .  has  the 
value 


S=‘jrpC2  :  -  (17) 

where  C  is  a  constant.  The  vortex  strength  is  given  by  the  pum 
of  the  tangential  velocities  on  the  two  sides.;  1  Thus  , 


and,  therefore. 


l/l  i  x  .'''ckp 

VyBx  dx 


(18) 


Jx=-1 


Carrying  out  the  indicated  calculations  in  appendix  D  leads  to 


C  a 


■f  ’V 


tot 


U'T0  +  +  (J0  "  ifl) 


J1  -  Yo  -  1(J0  *  h) 
;Tl  +  •%'  +  1(J0  -'’Ti) 
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where--  means  jQ(k)  and  so  forth  and  Yfl(k)  and  Y^(k)  are 
Bessel  functions  of  the  second  kind  and  of  the  zero  and  first  order. 
The  real  part  of  this’  equation  is  . 


C  /h  vf3o(y(X  +  Y)  cos-  cut-.  -  _Y(X  -  Y)  sin  cut 


’) 


where 


X 


(j0  -  Yj)2  *  (h  *  Y0)2  . 


(Jo  -  Ti)g  -  CJ1  *  Yo>g 

_(?o  -  Yx)2  +  (*i  +  Yo)2]2 


and  finally 


B  =  T^ph^Po2  (X  +  Y  cos  2cut  -  Z  sin  Scot) 

1  ■  v  •!  .  “  '  ’ 

where 


(19) 
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UXNAMICS  OF  AN  AIRFOIL  IN  AN  OSCILLATING  STREAM 


The  problem  treated  herein  is  that  of  an  airfoil  executing 
torsional  vibrations  in  an  oscillating  stream.  Because  of  assump¬ 
tions  (c)  and  (d)  of  the  Introduction  -  namely,  that  the  oscillations 
are  small  and  that  the  wake  is  plane  -  the  aerodynamic  moments  due 
to  partial  motions  of  the  fluid  and  those  due  if®  motions  of  the 
airfoil'  can  be  treated  independently. 

The  pitching  moment  acting  on  an  airfoil  undergoing  angular 
oscillations  in  a  uniform  stream  is  (reference  5) 


Thus,  the  total  aerodynamic  moment  acting  on  an  oscillating 
airfoil  in  an  oscillating  stream  is  the  sum  of  equations  (15) 
and  (20).  The  equation  of  motion  is  obtained  by  expressing  the 
equilibrium  of  the  aerodynamic  moment,  the  moment  of  inertia, 
and  the  mechanical  restoring  moment.  Thus,  if  structural  damping 
is  neglected,  the  equation  of  motion  Is 


Aerodynamic  moment  »  Iad  +  Raa 

where 

Ia  moment  of  Inertia  about  x  =>  a 
Ra  . torsional  stiffness  .  - 

Combining  equations  (15),  (20),  and  (21)  gives 


£a2  +  +  a2)J  &  +  1-2  (a  +  |)c(k) 

+  Op  ■  $  yS  (a  +  a)cW]“ 


a 


_  2  KVe 

b2 


(a  +  |c(k)(j0  -  iJ-i)  +  uj< 


icnt 


(21) 


(22) 
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m 


& 

-P 

I 

© 

© 

a* 

■ft 


OJ  © 

cm  .a 

-p 

d  <h 
o  o 

I  -H 

-P  © 


CM 


■sea1 


CVJ 
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where  F(Jc)  +  iG(h)  =  C(k) ,  TheoAoreen  s  C -function. 
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DISCUSSION  OF  RESULTS 


Propulsive  force.-  It  has  "been  shown  (see  equation  ( 19) )  that 
the  propulsive  force  acting  on  a  fixed  airfoil  at  zero  anglp  of 
attack  in  an  oscillating  stream  is 


■s  a  itphv^Q^(X'+  Y  cos  2u>t  -  Z  sin  2pt) 


Since  2ifpbv%Q  is  the  .stationary  value  of  the  lift  on  an  airfoil 
at  angle  of  attack  pQ,  let 


Lq  «=  2rtphv2 


Then 


B  si 


*oPo 


(X  f  Y  cos  Opt  -  Z  sin  2tot) 


(25) 


and  the  average  value  of  s  is 


LnP, 


r0p0 


In  figures  3  and  4  curves  are  presented  that  show  the  variation 
with  wave  lengths  of  stream  oscillations  of  the  coefficients  X,  Y, 
and  Z  appearing  in  the  equation  for  the  propulsive  force.  For  very 
lew  frequencies  -  that  is,  for  long  wave  lengths  of  stream  oscilla¬ 
tions  -  X  and  Y  approach  ’the. value  1,  and  Z  "becomes  zero.  Thus, 
as  k—?0,  equation  (25)  "becomes 


■SR# 


o  0 
2 


(1  +  cos  2 


(2 6) 


This  resitLt  is  exactly  that  which  is  to  "be  expected  from  quasi- 
statlonary  considerations  in  which  the  lift  is  assumed  to  "be 
instantaneously  that  value  prescribed  by  the  geometrical  angle  of 
attack^  that  is,  the  shed  vorticity  produced  by  variations  in 
angle  of  attack  is  assumed  to  appear  instantaneously  at  a  point 
infinitely  distant  from  the  airfoil.  Thus  (see  fig.  5)  the 
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perpendicular  force  L  on  an  airfoil  at  angle  of  attack  $  with 
respect  to  the  stream  is 


L  =  2rtpv^b(3 


and  because  from  well-known  considerations  of  two-dimensional 
flows  there  can  be  no  induced  drag,  the  propulsive  force  must  be 


s  =  LP  «  2npbv%^ 


and  if  p  =  p0  cos  cot,  . 


a  a  2ffpbv20Q^cos‘T»t 

»  L()3o(|  +  |  C0B 


which  is  the  same  as  equation  (26) . 

The  reason  for  plotting  the  .coeff icients  Y  and  Z  against  k 
as  well  as  against  ^ ,  where  is  the  wave  length  of  stream 


oscillations,  is  that  in  the.  neighborhood  of  * 

*  ’  Jc 


0  the  values  of 


Y  and  Z  fluctuate  infinitely  many  times.  This  behavior  is 
brought'  out  in  figure  4  where  it  may  be  seen  that  between  . 


k  *»  1 


t-1 


and  k  =03  (i  s  o)  the  curves  oscillate  about 


zero  with  decreasing  amplitude. 


Forced  vibrations.-  An  analysis  is  made  of  the  response  (see 
equation  - ( 2b))  of  an  airfoil  in  an  oscillating  stream  with  particular 
emphasis  on  the  parameters  a  and  Ra,  these  parameters  having 
qualitative  as  well  as  quantitative  effects  on  the  values  of  the 
response.  .The  stiffness  in  torsion  Ba  is  related  to  the  natural 
frequency  V'  '  for  zero  stream  velocity  by  the  equation 
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where  I„ 

cl 

airfoil . 


+  npb*^  +  a2^  is  the 
From  equation  (27), 


total  moment  of  inertia  of  the 


B 


a 


2ftmv2 


+ 


where  k'  = 

v 


In  figure  6  are  shown  response  curves  for  an  airfoil  hinged  at 
the  leading  edge  (a  =  -1.0.)  and  having  various  natural  frequencies 
of  oscillation.  The  following  values  of  the  parameters  k  and 
were  chosen  as  being  within  the  practical  range  of  application: 
k  =  0.0653^  ra2  =  i  (flat- plate  mass  distribution).  The  response 

for  any  natural  frequency  is  never  very  large  (at  most  about  —  n  2.5)* 

\  Pp  ' 

and  the  response  decreases  with  Increasing  stiffness.  The 
value  of  the  stream  frequency  at  which  maximum  response  occurs  is 
seen  to  correspond  more  closely  to  V 1  as  the  natural  frequency 
increases.  Even  for  a  freely  hinged  airfoil  (k1  p)  '  a  sort  of  . 
resonance  frequency  exists.  (See  fig.  J.)  Thus  the  use  of  such 
a  device  for  measuring  angular  variations  in  an  oscillating  stream 
would  be  valid  only  in  the  range  of  very  low  reduced  frequency  (long 
wave  lengths)  in  which  the  response  approaches  unity. 

Somewhat  different  phenomena  occur  when  the  hinge  point  is 
behind  the  quarter- chord  point  of  the  airfoil.  In  general  it  may 
be  stated  that  the  response  is  greater.  In  particular.^  critical 
stiffness  exists  below  which  the  airfoil  is  in  unstable  equilibrium. 
The  condition  for  divergent  motions  of  the  airfoil  is  that  the 
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IT 


coefficient  of  a  in  equation  (22)  "be  less  than  or  equal  to  0 
for  k  *  0.  The  critical  stiffness,  therefore,  is  defined  by 


B. 


2Kmv2 


=  a 


♦i 


In  figure  8  are  presented  some  results  for  the  hinge  placed 
near  the  center  of  the  airfoil.  The  values  of  the  parameters  a, 
k f  and  i*a2  chosen  were:  a  =  -0.1,  «  =  O.OS^S,  and  ra^  = 

The  reduced  critical  frequency  for  this  case  is  k’  =  0 .385 .  For 
stiffness  values  somewhat  higher  than  the  critical,  the  response 
is  not  unduly  large  and  again,  as  when  the  hinge  was  at  the 
leading  edge,  the  maximum  response  with  large  values  of  the 
stiffness  occurs  at  a  stream  frequency  close  to  the  frequency  v 1 , 


CONCLUSIONS 


The  velocity  potential,  lift  force,  moment,  and  propulsive 
force  on  a  two-dimensional  airfoil  in  a  stream  of  periodically 
varying  angle  of  attack  has  "been  derived  on  the  "basis  of  non¬ 
stationary  incompressible  potential-flow  theory  which  includes 
the  effect  of  the  continuous  sheet  of  vortices  shed  from  the 
trailing  edge.  Application  of  these  results  was  made  in  an 
analysis  of  the  variation  with  frequency  of  the  propulsive  forces 
on  an  airfoil  in  an  oscillating  stream  and  in  an  analysis  of 
the  problem  of  forced,  vibration  of  an  airfoil  in  an  oscillating 
stream  with  consideration  of  the  stiffness  of  the  airfoil  and 
the  position  of  its  torsion  axis.  The  following  conclusions 
were  indicated: 

1.  The  value  of  the  propulsive  force  acting  on  an  airfoil 

in  an  oscillating  stream  is  sufficiently  large  to  be  of  practical 
importance . 

2.  The  amplitude  of  vibration  of  an  airfoil  in  an  oscillating 
stream  is  critically  dependent  on  the  stiffness  of  the  airfoil 
and  the  position  of  its  torsion  axis.  In  general,  amplitudes  of 
vibration  are  smaller  when  the  torsion  axis  is  ahead  of  the 
quarter-chord  point  and  larger  when  the  torsion  axis  is  behind 
the  quarter -chord  point.  Because  of  the  nonuniform  response  of 
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a  freely  hinged  airfoil  the  use  of  such  a, device  for  the  measurement 
of  angular  variations  in  an  oscillating  stream  would  "be  restricted 
to  the  range  of  very  low  frequency  in  which  the  response  approaches 
unity., 


Langley  Memorial  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  April  28,  19^7 
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APPENDIX  A 

NONdRCDLATOEy  VELOCITY  POTENTIAL  FOE  AN 
AIRFOIL  IN  A  SINUSOIDAL  GUST 


The  problem  of  finding  the  velocity  potential  for  an  airfoil 
having  a  certain  normal  velocity  distribution  is  solved  by  the 
method  indicated  by  von  Karman  end  Burgers  (reference  6,  p.  44)  • 

Represent  the  wing  by  a  circle  (fig*  9).  Place  a  source  of 
strength  2e  at  the  point  ,  y-jJ  on  the  circle  and  a  Bink  of 
strength  -2e  at  -y^)  .  The  velocity  potential  of  this 

source- sink  pair  in  the  notation  of  the  present  paper  iB  given  by 
(reference  5) 


£p  = 


6 

2  it 


(x  -  x^)2  +  (y  -  yp)2 
log  - - - - - - - - - 

<*  -  xi)2  +-(y  +  xi)2 


The  transformation  of  the  circle  to  its  diameter  is 


y  =  t/l  -  x2,  x  =  x 
For  the  airfoil  in  a  sinusoidal  gust 


€  =  vfJoeiC®*'**!) 


Thus 


£cp  =  ZES.  eicjte-ikx! 
2ir 


log  -  *i)2  +  (y  -  yi)2 


(x  -  xx)2  +  (y  +  yi)2 
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Integration  "by  parts  leads  to 
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But,  from  relation.  (4)  on  paqe  48  of  reference  7, 


•where  Jq(1e)  is  a  Bessel  function  of  the  first  kind  and  zero  order. 
Therefore, 


|£  *  i*J0(k)  -  ixf  (A2) 

Equation  (A2)  is  a  non -homo geneoue  differential  equation  of'  the  first 
order.  The  homogeneous  part 


8f 

St 


-ixf 


has  as  solution 


f  » 


ce 


-ikx 


where  c  is  an  arbitrary  constant  with  respect  to  k.  The 
particular  solution  is  obtained  by  the  method  of  variation  of 
parameters  (reference  8,  p.  114) ,  Let 


f  =  g(k,  x)e"ikx 


(A3) 


-ikx 


Then 


8f  =  e-ikx  -  ixge 
8k  8k  ^ 
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Combining  this  expression  with  equation  (A2)  gives 

fjj  - 

Integration  of  which  leads  to 


g  =  i it  f  eiuxJo(u)  du 


Combining  equations  (A3)  and  (Ah)  gives 


f(k^x)  =>  irte“ik3:  eiuxJ0(u)  du 


which,  when  substituted  into  equation  (Al),  gives  for  the  non- 
circulatory  velocity  potential 


(A4) 


NACA  TN  No.  1372 


r\  <*- 


APPENDIX  B 

LIFT  AND  MOMENT  FOB  NONCIRCULATOBY  FLOW 
Lift 


The  total  force  on  the  airfoil  is  given  "by 


p'b£ 


Ap  dx 


vhere,  from  equation  (2) 


Ap  -  -2p  fz  ^ 

bx  at 


But 


Therefore 


**-1  =  <P*--1  *  0 


P  =  -2ph 


f  1  & 

J-i  at 


dx 


From  equation  (l). 


I  *  -  ^Tr2.  A  -  X2  e^-itac  /  eltncJ0(u) 


du 


(Bl) 


Thus 


P 


iSph^vg^o 


.loot 


k 


/^  0-  -  x2  e"iloc  dx/  eiuxJ0(u)  du  (B2) 


2k 
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Let 


fL(k)  o  J  /l  -  x2  eTitac  Ax  J*  eluxJ0(u)  du  (B3) 


Interchanging  the  order  of  integration  gives 


fx(k)  «  J  J0<u)  du  iy 


2  i(u~k)x 


x  e 


4x 


From  reference  7,  page  48, 


*(.)  -  -  AL <  <  f1  a  -  t=)v-l  6lztdt 


Therefore 


(B4) 


(B5> 


J  (U  -  k)  «  ■  liLl  *1  f  1  (1  -  X2)l  0i(u-k)x  dx 

»  f  1  (1  -  *2)1  di  • 

TT 


Substituting  this  expression  into  equation  (B4)  gives 


fl(k)  J{ C  J°(u)  JiCu " k> 
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Then 


u  -  k  =  -w 


du  «  -dw 


^(k)  -  *  j  ^(-v)  j0(r  - »)  t 

k 


=  -n  f  :J- -  w)  52  . 

WQ  •*•  u  w 


*£  Vv>' 


J-0(fc  -  v)  % 


But  from  reference  7,  page  380/ for  u  >  0  and  v.>  -1 


S  j„(t)  t)  £s.5ai 

J0  ^  t  u 


Therefore 


fx(k)  =  «Jx(k) 
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and  equation  (B2)  ’becomes 


^(k)  .tot 


(B8) 


Moment 

The  noncirculatory  moment  about  the  point  x  =  a  (fig.  2) 
is  obtained  from  the  integral 


Ma  «  b£ 


J'  Ap(x  -  a)  dx 
which,  combined  with  equation  (2),  gives 


Ma  a  2pbv  J  cp  dx  +  2pb2a  J  dx  -  2pb2  J 


By  use  of  equations  (l)  and  (B3.),  this  equation  becomes 


2pb2v2pn  nl  r - -  r* 

=  -  — _2  etot  l/l  -  x2  e’ikx  dx  £  elux  JQ(u)  du 


„iux 


-  g^Wgp  _itut  r  1 


k 


ft ~ 


x&  e"iJa:  dx  j  elux  J0(u)  du  (B9) 


/ 

J0 


iux 


2ipb-WpQ 


,in)t 


r 


x2  e-ikx  4x  /  eiux  j  (u) 
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The  first  two  integrals  on  the  right  side  -nf  this  expression  are 
already  known  (see  equations  (B3)  and  (B7) )  •  In  order  to  obtain 
the  third  integral,  let 


fo(k)  »  f  1  \jl  -  x2  e"1^  dx  f  e*113-  Jo(u)  (BIO) 

^1  J0 


Interchange  the  order  of  integration.  Then 


f2(k) 


J0(u)  du 


Qi(u-k)x 


dx 


(Bll) 


But 


ei(u-k)x  ^ 


|'(1  -  x2)3/2  ei(u-k)x 


.1 

-1 


(B12) 


+ r  u  -  z2)3/2.1^*  ax 

"3  *1 

*  ^L—^l  f  (1  .  3c2)3/2  dx- 

3  -  -  -1 


and,  from  equation  (B5) 


-  x2)3/2  dx  = 

(u  -  -k)2 


(B13) 


28 


MCA  TN  No*  1372 


Combine  equations  (B12)  and  (B13)  and  substitute  into  equation  (Bll) . 
Then 


fo(k)  »  i«  f.  J0(u)  J2(u  -  k)  — ~ 

G  i'  ft  .  U  ~  K 


Let 


and 


Then 


u  -  k  =*  -v 


du  =  -dw 


f2(k)  ®  in  jT  Jg(-w)  JQ(k  -  w)  «*  -iif  J2(-w)  JQ(k  -  w)  ~ 


-,-i*  Jg(w)  J0(k  -  w)  ^ 


which  becomes  (see  equation  (B6)) 


f„(k)  =  -  is  Jp(k) 


(BlU) 


The  equation  for  the  moment  about  x  «  a  can  now  be  written  as 


**a 


2Ttpb2v23n  .  .  2itipb3aovf50  ,  ^ 

- - - £  J3_(k)  etot - r - 2  J  (*)•*«* 


k 

7tpb3a)vp0 


(B15) 


J2(k)  etot 
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The  recurrence  formula  from  reference  7 ,  page  1? , 


Jn-l(k)  +  ^n+l^  =  ^  ^n(^) 

gives  for  n  =  1 

J2(tf  -  -Jo(k)  +  |  Jl(k) 


Substituting  this  expression  into  equation  (B15)  and  making  use  of 
the  definition  k  =  ~~  yields 


Ma  =  -2nipb2av2pQ  ^(k)  eijjt  •. -itpb^Po J0(  k)  'etot  '  (Bl6) 


APPENDIX  C 


EVALUATION . OP  Q  (EQUATION  ( 8) ) 


Equation  (8)  is 


Q  =  -  2°  ei(cbt-k)  eiu  j  (u)  du 

k  '■'o  -  '~ 

^0  lcot  fk  i(u-k) 
s*  -  - .  e  /  ©  Jn(u)  du 

k  J0 


But 


oos  (u  -  k)  J0(u)  du 


+  i 


sin  (u  -  k)  J0(u)  du 


From  reference  7>  pageB  380  and  381, 


J k  cos  (k  -  u)  Jq(u)  du  =  kJQ(k) 


J'  sin  (k  -  u)  J0(u)  du  =  kj^(k) 
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Therefor© 


f  el(u'k) 


J0(u)  du  =  kT0(k)  -  ikJiCk) 


=  k[j0(k)  -  U^k)] 


and 
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APPENDIX  D 

CALCULATION  OP  THE  PROPULSIVE  FORCE 
From  equation  (l) 


^0  ei(oit+k) 
k 


e‘lu  J  (u)  4u 


e'1(u‘S)  J0C»)  to 


The  integral  on  the  right  aide  of  this  equation  is  the  complex 
o on Jugate  of  that  in  appendix  C.  Thus 


J0(u)  du  =  k(JQ  -  ij^) 
“  KJo  +  lj,i) 

where  Jq  means  jQ(k),  and  so  forth.  Therefore 


(Mh-i m  ■fe4?5,.l  ^ (Jo  +  1Jx) 


Now,  for  the  velocity  due  to  the  wake, 


.ioit 


(ML) 
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"but  (see  equations  (7)  and  (8)) 


JL_ 

2jt 


Xq2-  1 


TJ  dZQ  «  Q, 


Therefore 


From  reference  7,  page  180, 


Bin  ^k:  coBh  t  -  i  n  ir^  cosh  nt  dt 


and 


Y  (k)  =  -  -  /  cos  (k  cosh  t 
n 


^k  cosh  t  -  |  n  it^ 


cosh  nt  dt 


(D2) 
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which  lead  to,  when  cosh  t  Is  replaced  by  Xq, 

^.5  r 21*2*0 

J 1  /*C)S  -  1 


Jl(k) 


XqCos  k'^Q  dXQ 


Y0U)  * 


cos  k  Xq  Ixq 
l/xQ2  -  1 


Equation  (D2)  "becomes  now 


XqsIh  k  xq  dXQ 

V'xq2  -  1 


(fi  '  Yo)  -  *(*>  *  Yi) 
1/2  CJi  +  Yo)  +  i(Jo  -  Yi) 


and  after  substitution  of  the  value  of  Q  from  equation  (C2) 
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t~)  =  -  (-  •1--^  (J0  -  iJi) 

vpz  / x=-l  (i/l  +  x/x=-l  i/2 


(ft  -  Yo)  -  KJo  +  YJ 


.loot 


LCJ1  +  Y0>  +  i(J0  -  Yi)j 
/ ••  ■  :  -  (D3) 

Combining  equations  (Dl)  :an&  (D3)  withr'(see  equation  (18)) 


C  = 


k=-l 


gives 


‘/I  Tf30etot  (J0  +  Ul)  +  (Jo  "  Ul)  ^J°  + 

I  *-  (  T_  j.  V-  ^  a.  A  (  _  V,  \ 


r  „ 


(Jl  +  ?0.)  +  HJ0  "  Yl) 

W 


=  -i/^  vp0(co8  cut  +  i  sin  <Dt)<(Jo 


1  + 


Jl2  -  jo2  +  yi2-yo2 


i- 


2Jl(J0Jr  +  YpYx) 

(Jl  +  Y0<)2  +  (J0  -  Tl)J 


+  iJx 


(ft  +  Y0)2  +  (J0  -  Y^ 


Jlg  -  Jp2  +  ft2  -  Y0g 

(Jl  +  Yo)g  +  (Jo  -  ft)2 


_i  2Jq(JqJi  +  YpYj) 

(ft  +  Yo)g  +  (Jo  -  ft)2 

the  real  part  of  which  is 


C  »  \fb  vp 


0 


i/^X_+_Y)  cos  cot  -  \f(T  -  Y)  sin  cot 


m) 
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where 


Wi  -  m)% 

(Jo  -  Tl)2  *  (Jl 


■  (£.f _ _ 1 - 

W  (Jo  -  Yl)s  +  (Jl'+ 

and 

Y  n  CJ0*1  ~  J^g  &J0  •  Yl.)g  -  *  h)fl 

f(Jo  -  Yi)a  +  (J1  -  Vf 

.  (S.f  (Jp  -  *lf  -  (Jj  ♦ 

vtk/  r  1 2 

.  [(J0  -  Y^2  +  (Jj.  +  Y0)2J  • 

and  use  has  "been  made  of  the  formula  (reference  7 ,  P»  77) 


J0T1  •  J1Y0  -  -  i 


Substituting  equation  (D5)  into  (see  equation  (17)) 


s  =  JtpC® 
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gives 


s  =  irpl>v2302 


(X  +  Y)  cos^ot  +  (X 


-2 


sin  cot  cos  cut 


,2.  2 


=  itpbv^Po*  (X  +  Y  cos  .2  cut  -  l/x 


Letting  Z  =  |/x2  -  Y2  .  gives  finally 


s 


jrptv^o2  (X  +  Y  cos  2  cot  - 


-  Y)  sin^rot 


-  Y2  sin  2  cot 


Z  sin  2  cot) 
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Fig.  3 
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f~/'gure  -4  >  ~~  Variation  with  reduced  fregucnc  y  of  efream  oscii/af/ons 
of  the  coe fficien is  Y  and  if  . 
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Fig.  5 
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ffgc/re  7 .  —  \rt/r/af/'on  rr/'/h  war  e  /ength  of  stream 
oae/tfah/ons  of  the  response  of  an  a/rfo/f  free/ g 
h/nged  at  }fs  /eac//r?g  edge. 
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d/gure  Q.  —  Is&r/at/on  with  rec/ucec/  frequence/  of 
-stream  osc/  i/a  f tons  of  the  response  of  an  a/rfo/Z 
hinged  near  its  center. 
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